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ABSTRACT 

We present a systematic study of the spectral consequences of departures from chemical equilibrium 
in the atmospheres of L and T dwarfs, and for even cooler dwarfs. The temperature/pressure profiles 
of the non-equilibrium models are fully consistent with the non-equilibrium chemistry. Our grid 
of non-equilibrium models includes spectra for effective temperatures from 200 K to 1800 K, three 
surface gravities, four possible values of the coefficient of eddy diffusion in the radiative zone, and three 
different CO/CH4 chemical reaction prescriptions. We also provide clear and cloudy model variants. 
We find, in keeping with previous studies, that there are essentially only two spectral regions where 
the effects of departures from chemical equilibrium can influence the predicted spectrum. These are in 
the M (^4—5 /im) and N (8—14 /im) bands due to CO and NH3, respectively. The overabundance of 
CO translates into flux suppressions of at most ~40% between effective temperatures of 600 and 1800 
K. The effect is largest around T c g rs 1100 K. The underabundance of ammonia translates into flux 
enhancements of no more than ^20% for the T e g range from 300 to 1800 K, with the largest effects 
at the lowest values of T ff. The magnitude of the departure from chemical equilibrium increases 
with decreasing gravity, with increasing eddy diffusion coefficient, and with decreasing speed of the 
CO/CH4 reaction. Though these effects are modest, they lead to better fits with the measured T 
dwarf spectra. Furthermore, the suppression in the M band due to non-equilibrium enhancements in 
the CO abundance disappears below ~500 K, and is only partial above ~500 K, preserving the M 
band flux as a useful diagnostic of cool atmospheres and maintaining its importance for searches for 
the cooler brown dwarfs beyond the Ts. 

Subject headings: stars: abundances — stars: atmospheres — stars: individual (Gliese 570D) — stars: 
low-mass, brown dwrafs 



1. INTRODUCTION 

It has been long realized that departures from the local 
chemical equilibrium (LCE) may play a significant role 
in the atmospheres of the substellar-mass objects (SMO), 
that is of giant planets and brown dwarfs. Essentially, 
some chemical reactions occurring in these atmospheres 
may be very slow, so that vertical transport via con- 
vective motions or eddy diffusion can lead to departures 
from LCE. The mechanism was first suggested to oper- 
ate in the Jovian planets of the solar system (Prinn & 
Barshay 1977, Barshay & Lewis 1978, Fegley & Prinn 
1985, Yung et al. 1988, and Fegley & Lodders 1994). 

Fegley & Lodders (1996) suggested that the same 
mechanism can also operate in the atmospheres of brown 
dwarfs. Noll, Geballe, & Marley (1997) studied in more 
detail the non-equilibrium enhancement of CO in the 
atmosphere of Gl 229B. Subsequently, Griffith & Yelle 
(1999) and Saumon et al. (2000) provided a more de- 
tailed analysis. The latter paper also suggested that 
a non-equilibrium depletion of NH3 may be important. 
Saumon et al. (2003) demonstrated the effects of non- 
equilibrium chemistry for a range of effective tempera- 
tures (T e ff between 800 and 1600 K). Recently, Saumon et 
al. (2006, 2007) showed that the non-equilibrium effects 
in the carbon and nitrogen chemistry have significant ef- 
fects on the predicted spectrum of brown dwarfs, and in 
turn fit the observed spectrum of the T7.5 dwarf Gliese 
570D (Saumon et al. 2006) and 2MASS 0415 (Saumon et 
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al. 2007) much better than models based on local chem- 
ical equilibrium. Finally, Leggett et al. (2007) extended 
the Saumon et al. (2006) models to a larger range of 
parameters, and presented a comparison of model pre- 
dictions with observed infrared colors. 

In this paper, we extend and generalize the Saumon et 
al. and Leggett et al. treatments by computing a large 
grid of solar-metallicity models for effective temperatures 
T c ff covering the whole domain of T and L dwarfs (700 
- 1800 K), as well as for cooler dwarfs (200 - 700 K), 
for three values of surface gravity (log g = 4.5, 5.0, 5.5 
cm s -2 ), for several values of the coefficient of eddy dif- 
fusion (see §3.1), and for three different reaction rate 
timescales for the carbon reactions (see §3.1). We com- 
pute self-consistent models in the sense that the atmo- 
spheric structure (T/P profile) is computed taking into 
account departures from chemical equilibrium. We have 
computed two grids of models, one cloudless, and one 
including clouds composed of silicate condensates. 

In §SJ we outline our procedures for computing model 
atmospheres, while allowing for departures from chemical 
equilibrium. In 53 we present and analyze the predicted 
spectra computed with and without chemical equilibrium 
for L and T dwarfs. We explore the dependence on T e g, 
gravity, and eddy mixing coefficient in the radiative zone 
and compare clear and representative cloudy models. We 
also discuss the effect of non-equilibrium chemistry for 
objects cooler than the coolest T dwarfs down to effective 
temperatures of 200 K. Then, in §31 we briefly present a 
comparison between theory and observation in the mid- 
IR for Gliese 570D. Finally, in $5j we summarize our 
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conclusions. 

2. MODELING PROCEDURES 

2.1. Non- equilibrium Carbon and Nitrogen Chemistry 

The chemistry of carbon and nitrogen in the atmo- 
spheres of SMOs is essentially described by the net reac- 
tions 

CO + 3H 2 < — >CH 4 + H 2 0, (1) 



for carbon, and 



N 2 + 3H 2 



2NH 3 , 



(2) 



for nitrogen. Because of the strong C=0 and N=N 
bonds, the reactions ([I]) and ((2]) proceed much faster from 
right to left than from left to right. For instance, for car- 
bon the reaction in which CO is converted to CH4 is very 
slow, and, therefore, CO can be vertically transported in 
the atmosphere to the upper and cooler regions by con- 
vective motions or eddy diffusion. The net result is an 
overabundance of CO and N 2 and an underabundance 
of NH3 and CH4 in the upper atmosphere. Since more 
oxygen atoms are now tied in CO, the mixing ratio of 
water is also reduced because of conservation of the total 
number of oxygen atoms. 

The mixing (vertical transport) timescale, i m i X) was 
suggested by Griffith & Yelle (1999) to be parameterized 

by 

W = H 2 /K zz , (3) 

where H is the pressure scale height, and K zz is the 
coefficient of eddy diffusion. Equation ((3|) applies in the 
radiative (convectively stable) zone. In the convection 
zone, the mixing time is given by (Saumon et al. 2006): 



^mix — 



3-ffc/Vc 



(4) 



where H c is the convective mixing length (adopted here 
to be equal to H) , and v c is the convective velocity. We 
stress that while the mixing time is well defined in the 
convection zone, it can only be roughly parameterized 
in the radiative zone. As discussed e.g. by Saumon et 
al. (2006, 2007), the reasonable values of K zz are in the 
range 10 2 - 10 cm 2 s -1 . 

The chemical reaction timescale, < c hom, generally de- 
creases rapidly with increasing pressure and tempera- 
ture. In the atmospheric regions where t m i x > i c hemj 
the species acquire their equilibrium mixing ratios. This 
occurs at temperatures and pressures larger than that 
at the depth where t mlx = t c \ lcm . At smaller pressures, 
the species acquire mixing ratios equal to those at the 
t m ix = tchem layer, as first suggested for the case of the 
CO/CH4 reaction by Prinn & Barshay (1977). 

The non-equilibrium mixing ratios thus depend sensi- 
tively on the chemical reaction timescales, which unfor- 
tunately are poorly known. For the N 2 — > NH3 conver- 
sion, we adopt the timescale given by Lodders & Fegley 
(2002): 



^chcm — 



1 



k N2 N(H 2 ) 

>vn, s.W x 10~ 8 exp ( -^^\ , (5) 



T 



where kn 2 is the rate constant for the N 2 — * NH3 reac- 
tion, and 7V(H 2 ) is the number density of H 2 . 



The reaction timescales for the CO — > CH4 conver- 
sion are uncertain because the correct chemical pathway 
is still unknown. We will, therefore, use three different 
timescales, suggested earlier in the literature: 
i) the "slow" chemical timescale after Prinn & Barshay 
(1977): 



chcm 



tco 



JV(CO) 



kco = 2.3 x 10~ 10 exp 



K C oW(H 2 )iV(H 2 CO) 
36200^ 



T 



(6) 



where N(M) is the number density of species M. 
ii) and iii) the "fast" chemical timescale after Yung et 
al. (1988) or Dunning et al. (1984). In both cases the 
timescale is given by 



iV(CO) 



(7) 



k 7 N(R) N(R 2 CO) ' 

where the rate coefficient differs in each reference. We use 
tabulated values presented in Yung et al. (1988 - their 
Table III), and interpolate to current temperatures. In 
this paper, we call the Dunning et al. timescale "fasti," 
and the Yung et al. timescale "fast2." As we shall show 
in §3, these two fast timescales produce very similar re- 
sults. 

2.2. Atmosphere and Spectrum Modeling 

We use the updated code CoolTLUSTY, described 
in Sudarsky, Burrows, & Hubeny (2003), Hubeny, Bur- 
rows, & Sudarsky (2003), and Burrows, Sudarsky, & 
Hubeny (2006 — BSH), which is a variant of the univer- 
sal atmospheric code TLUSTY (Hubeny 1988; Hubeny 
& Lanz 1995). CoolTLUSTY solves self-consistently a 
set of radiative transfer equations for selected frequency 
points (typically 5000 points logarithmically spaced be- 
tween 7 x 10 14 Hz and 10 12 Hz), and the equation of 
radiative+convective equilibrium. The atmosphere is as- 
sumed to be in hydrostatic equilibrium. The opacities 
are not computed on the fly; instead, they are interpo- 
lated from pre-calculated opacity tables for the current 
values of temperature and density. Our opacity database 
is described in detail by Sharp & Burrows (2007). The 
atmospheric structure is computed iteratively, applying 
the Newton-Raphson method (first applied in the case 
of stellar atmospheres by Auer & Mihalas 1969 with the 
name "complete linearization" ) . The discretized and lin- 
earized equations are re-organized in the so-called Ry- 
bicki scheme (Mihalas 1978), which renders the com- 
puter time linearly proportional to the number of fre- 
quencies (and not proportional to the cube of the num- 
ber of frequencies as in the case of the original complete 
linearization). This allows us to avoid using the acceler- 
ated lambda iteration (ALI) scheme to treat the radiative 
transfer equation, which leads to faster and more stable 
convergence, without compromising on computer time. 

We have found that the most stable procedure to com- 
pute non-LCE models is to perform several iterations 
(typically 6) of the linearization method assuming chem- 
ical equilibrium (using the algorithm of Burrows & Sharp 
1999 to determine the equilibrium abundances of all 
species), and then to switch the non-LCE algorithm. In 
each subsequent iteration, one calculates the intersection 
point where the mixing time for the current T/P profile 
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equals the reaction time (computed again for the current 
T/P profile); the number densities of CO, CH4, H2O, 
and NH 3 are set to constant values equal to those found 
at the intersection points for pressures lower than the 
pressure at the intersection. The opacities in these re- 
gions are modified accordingly. The next iteration of 
the linearization scheme is then performed with the new 
opacities. The overall process usually converges in a few 
additional iterations. 

When clouds are taken into account, we employ the 
scheme described in BSH. In the case of L and T dwarfs, 
we represent individual cloud decks of many condenstate 
species by a single extended cloud deck of one representa- 
tive condenstate. As in BSH, we take forsterite as a rep- 
resentative condensate, and assume the cloud deck to ex- 
tend between the intersection of the condensation curve 
of forsterite with the current T/P profile, and the pres- 
sure where the local temperature equals 2300 K (which 
schematically represents the condensation curve of the 
most refractory condenstate), and with exponential de- 
creases of particle density on both sides of the cloud. In 
the terminology of BSH, this is an E-type cloud. In all 
the present models, we assume the modal particle size 
100 fim. 

3. RESULTS 

3.1. L and T Dwarfs 

Our effort has generated an extensive grid of non- 
equilibrium model spectra for effective temperatures 
from 700 K to 1800 K (in steps of 100 K), for three values 
of surface gravity, g = 10 45 , 10 5 , and 10 5 5 cm s~ 2 , for 
four values of the coefficient of eddy diffusion in the radia- 
tive zone K zz = 10 2 , 10 4 , 10 6 , and 10 8 cm 2 s~\ and for 
three different CO/CH4 reaction time prescriptions 2 . We 
have also provided clear and cloudy model variants. The 
latter assumes an E-type forsterite cloud with a modal 
particle size of 100 /im. In the present paper, we consider 
only solar-composition models, where the solar elemen- 
tal abundances are from Asplund, Grevesse, & Sauval 
(2006). 

We will first discuss composition profiles and their 
trends with effective temperature, surface gravity, diffu- 
sion coefficient, and the CO/CH4 reaction speed. Figure 
[T]displays the composition profiles of the non-equilibrium 
species, CO, CH4, NH 3 , and H 2 0, as a function of the 
local temperature in the atmosphere, for a representa- 
tive model with T cff = 900 K, g = 10 5 5 cm s" 2 , and 
K zz — 10 4 cm 2 s _1 . The full lines show the non- 
equilibrium mole fractions or, equivalently, the number 
fractions (expressed as a fraction of the total number of 
particles), while the dashed lines show the equilibrium 
number fractions. For this model, the non-equilibrium 
effects in water and methane are negligible. The non- 
equilibrium number fraction of ammonia is smaller than 
the equilibrium one for T < 1200 K, as first shown by 
Saumon et al. (2006). The non-equilibrium mixing ra- 
tio of CO is significantly larger than the equilibrium 
one, which in turn leads to a significant strengthening 
of the 4.7-/im feature (and, to a lesser extent, to other 
CO features at 3.4 and 2.8 /im) in non-equilibrium mod- 
els. Moreover, the non-equilibrium mixing ratio of CO 

2 Electronic versions of these spectral models can be found online 
at http : //zenith. as . arizona. edu/ "burrows. 



is larger for the slower CO/CH4 reaction, because the 
intersection of the t m ; x (T) curve with the tco(T) curve 
occurs deeper in the atmosphere for higher temperature, 
at which the equilibrium CO mixing ratio is higher (see 
also Fig. [5J. Figure [T] is analogous to Fig. 3 of Saumon 
et al. (2006). 

To depict the trend of composition profiles with ef- 
fective temperature, we display in Fig. [2] composition 
profiles as in Fig. CD but for T eff /100 K = 8, 10, 12, 14, 16, 
and 18. To avoid cluttering the figure, the models are not 
labeled by the values of T e ff. Instead, the models are dis- 
tinguished by the fact that the higher the effective tem- 
perature is, the longer the curves continue toward high 
temperatures (i.e., to the right side of the plot). First, 
we stress that the equilibrium abundance patterns follow 
from the basic chemistry (Burrows & Sharp 1999; Lod- 
ders & Fegley 2002), and are well understood. Briefly, 
the same local temperature is reached at lower pressure 
for models with higher effective temperatures (see Fig. 
[7] for an explicit demonstration). Due to Le Chatelier's 
principle, at a given local temperature the abundance 
of CH4 decreases and the abundance of CO increases 
with increasing T e g. Analogously, the abundance of NH3 
decreases, and the abundance of N2 (not shown here) 
increases with increasing T c g. These effects are clearly 
shown in Fig. O 

For all models displayed here, the non-equilibrium ef- 
fects in water and methane are negligible. The non- 
equilibrium abundance of CO is always larger than the 

equilibrium abundance for T < 1000 K. For very low 
T g, although the non-equilibrium abundance of CO is 
significantly larger than the equilibrium abundance, it is 
still small in absolute value, and, therefore, does not have 
much of an effect on the model structure and predicted 
spectra. For the highest T e g, the equilibrium abundance 
profile of CO becomes flatter, and, consequently, the 
differences between the non-equilibrium and equilibrium 
abundance profiles of CO decrease, and so do the differ- 
ences in the predicted spectra. This is demonstrated in 
Fig. [5J For ammonia, departures from equilibrium be- 
come smaller for higher effective temperatures, because 
the equilibrium ammonia composition profile becomes 
flat. Therefore, fixing the ammonia number fraction at 
the intersection of the t mlx and the £nh 3 curves does not 
lead to a significant difference in the composition. There- 
fore, the non-equilibrium effects in ammonia, and, con- 
sequently, in the predicted flux in the 8-14 /im region, 
decrease with increasing T c g. This also is demonstrated 
in Fig. [8] Moreover, because of its diminished number 
fraction, ammonia becomes a less important source of 
opacity as the temperature increases. 

In Fig. El we show the trend of the composition profiles 
with the surface gravity, for T c g = 900 K, K zz — 10 4 cm 2 
s _1 , and the fast CO/CH4 reaction. The intersection of 
the i m ; x and the too curves occurs at essentially the same 
temperature for all three surface gravities (see Fig. |4|). 
However, at a given temperature, the model with the 
lowest gravity exhibits the lowest pressure, and, there- 
fore, the CO mixing ratio at the intersection is higher 
for the lower-gravity models. Consequently, the effects 
of departures from chemical equilibrium are greater for 
the low-gravity models. 

Figure |4] displays the mixing time and the reaction 
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times for the models displayed in Fig. [3l Notice the 
abrupt change of the mixing time at T « 1800 K; this 
corresponds to the onset of convection, where the mix- 
ing times are about five orders of magnitude faster. To 
explain the behavior of models with different coefficients 
of eddy diffusion, and with different reaction timescales 
for the CO/CH4 reaction, we display in Fig. [5] the mix- 
ing times for four different diffusion coefficients, K zz = 
10 2 ,10 4 ,10 6 , and 10 8 cm 2 s _1 . By increasing the coef- 
ficient of eddy diffusion, K zz , the mixing time decreases 
and, consequently, the intersection of the chemical reac- 
tion timescale curves and the mixing time curve occurs 
at higher temperatures. Therefore, the non-equilibrium 
mixing ratio of CO is higher. For the N2/NH3 reaction, 
and for K zz > 10 4 cm 2 s _1 , the chemical time curve in- 
tersects the mixing time curve in the convection zone at 
very high temperature (T « 2100 K). 

In Fig. [6l we display trends of the chemical reaction 
time with the effective temperature. The N2/NH3 reac- 
tion timescale is relatively insensitive to the model ef- 
fective temperature, while for the CO/CH4 reaction, the 
intersection of the curves moves to higher temperatures 
and pressures for higher effective temperatures. It would 
seem at first sight that the departures from chemical 
equilibrium would increase with increasing T c g. How- 
ever, as we have shown in Fig. [2j the abundance profile 
of CO becomes flatter and, thus, the departures from 
chemical equilibrium decrease instead. 

In Fig. [7] we display the temperature/pressure pro- 
files for selected models of our grid. The differences in 
the local temperature between the equilibrium and non- 
equilibrium models are not large, but nevertheless appre- 
ciable. 

We now turn to emergent spectra. In Fig. [8j we dis- 
play the emergent spectra for all effective temperatures, 
and for one representative value of g = 10 5 5 cm s~ 2 , for 
K zz — 10 4 cm 2 s _1 , and for the "slow" reaction rate of 
the CO/CH4 reaction of Prinn & Barshay (1977). There 
are essentially only two wavelength regions where the 
departures from chemical equilibrium influence the pre- 
dicted spectrum, the region around 4.7 microns (due to 
non-equilibrium abundances of CO), and a wide region 
between 8 and 14 microns, due to non-equilibrium abun- 
dances of NH3. The CO feature has its peak for effective 
temperatures between 800 K and 1200 K; the effect gen- 
erally decreases for hotter models and essentially disap- 
pears at 1800 K. The non-equilibrium effect of ammonia 

already disappears for T c g ~ 1200 K. The behavior of 
the associated spectra is readily understood upon exam- 
ination of the composition profiles displayed in Figs. Hi- 
ll] (see also the related discussion). Figure [8] is a synoptic 
display of the overall model trends, but is not meant to 
communicate details of the predicted spectra. Therefore, 
we later present in Figs. [TTJ[12l and H~3l enlarged plots of 
the two non-equilibrium regions, the CO region between 
3 and 5 /jm, and the ammonia region between 8 and 14 
fiva. In these and subsequent figures that display the pre- 
dicted spectra, we smooth the the computed spectrum by 
a 10-point boxcar average, which corresponds to a reso- 
lution of about 1/100. The corresponding models with 
clouds are displayed in Fig. [9] and behave analogously, 
thus demonstrating that the presence of clouds does not 
lead to any qualitative differences in the magnitude of 
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departures from chemical equilibrium. 

Figure [10] shows the sensitivity of the predicted spec- 
trum to the surface gravity. The non-equilibrium ef- 
fects generally increase with decreasing gravity, as was 
explained in the discussion of Fig. [3] Figure [IT] presents 
a more detailed view of the CO non-equilibrium region 
between 3 and 5 /im, for one particular model with 
T cff = 900 K and g = 10 5 5 g cm- 2 . We display both 
the LCE model and four non-equilibrium models with 
K zz = 10 2 , 10 4 , 10 6 , and 10 8 cm 2 s _1 . All models are 
for the "slow" reaction rate of Prinn & Barshay (1977) 
for the CO/CH4 reaction. An analogous plot for the 
models with the "fast2" reaction of Yung et al. (1988) is 
displayed in Fig. [T2] The magnitude of departures from 
LCE increases significantly with increasing K zz . This is 
because the mixing time decreases with increasing K zz 
and, thus, the intersection of the t miyi and the tco curves 
occurs deeper in the atmosphere, where the mixing ratio 
of CO is larger. This was already demonstrated in Fig. [5] 
As Fig. [T2] shows, the departures from LCE are smaller 
for the faster CO/CH4 reaction rate, as demonstrated 
previously in Fig. [T] 

In Fig. I13[ we display a more detailed view of the am- 
monia region between 8 and 14 fim, for the same models 
as displayed in Fig. 1121 The departures from chemi- 
cal equilibrium are essentially the same as predicted by 
Saumon et al. (2006, 2007). Unlike the case of CO, the 
abundance of ammonia and, therefore, the flux in the re- 
gion between 8 and 14 fim, is relatively insensitive to the 
value of K zz . This is because the abundance profile of 
ammonia is rather flat (see Figs. [T]-[3]), and, therefore, 
the depth of intersection of the £n 2 an d the i m i x curves 
does not matter very much. The departures of ammo- 
nia from LCE are largely insensitive to the speed of the 
CO/CH4 reaction, as expected. We do not show this 
in the figure because the corresponding plots are essen- 
tially identical. The predicted departures of NH3 from 
LCE seem to be quite robust. 

To better display departures from chemical equilib- 
rium, we plot in Fig. [T3] the ratio of the integrated flux 
between 4.4 /im and 5 /xm of the non-equilibrium and 
equilibrium models, as a function of effective tempera- 
ture. The upper panel displays the dependence of the 
non-equilibrium ratio for g = 10 5 cm s -2 , for the "slow" 
CO/CH4 reaction, and the "fast2" CO/CH4 reaction. 
Each set of models is computed for three values of K zz . 
The curves are labelled by the value of log 10 K zz in cm 2 
s _1 . As discussed above, assuming the "slow" CO/CH4 
reaction increases the departures from LCE significantly, 
reaching about 30% at T cS = 1100 K (at K zz = 10 4 
cm 2 s -1 ), and decreasing toward higher temperatures. 
For the "slow2" reaction, the maximum effect is around 
Tcff = 1300 K, where it reaches about 10% (again, at 
K zz = 10 4 cm 2 s -1 ). The lower panel displays the de- 
pendence of the non-equilibrium ratio on surface gravity 
for two sets of models. In all models, K zz = 10 4 cm 2 s . 
The upper set of models (diamonds connected by dashed 
lines) are for the "fast2" CO/CH4 reaction, while the 
lower set of models (stars connected by dotted lines) are 
for the "slow" CO/CH4 reaction. The curves are labelled 
by the value of log 10 g in cm s~ 2 . As discussed above, 
the lower-gravity models exhibit larger departures from 
equilibrium for temperature below 1300 K. However, the 
trend is reversed for T c g > 1300 K, where the lower grav- 
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ity models are closer to equilibrium. 

3.2. Beyond the T Dwarfs 

In this section, we briefly explore the effects on spec- 
tra of departures from chemical equilibrium for brown 
dwarfs cooler than the coolest known T dwarfs (see also 
Burrows, Sudarsky, & Lunine 2003). To this end, we 
extend the grid of models described in W3.ll and com- 
pute models for T ofr = 200, 300, 400, 500, and 600 K. 
For the models below 500 K, we include the effect of a 
water cloud with a modal particle size of 100 /im and a 
super-saturation parameter of 1.0%. The other basic pa- 
rameters are the same as in ^3.11 In Fig. [TBI we display 
the mixing time and the reaction times for three models 
with T c ff = 200, 400, and 600 K, for g = 10 5 cm s~ 2 , 
K zz = 10 4 cm 2 s _1 , and the "fast2" CO/CH 4 reaction. 
In all models, the intersection of the mixing time and 
the reaction times occurs in the convection zone. There- 
fore, the actual value of the diffusion coefficient K zz , 
which influences the mixing time only at temperatures 
and pressures lower than those at the intersection, is in- 
consequential. (In other words, K zz influences the T mrx 
curve only on the left of the intersection point, which 
is already determined by the local value of the convec- 
tive velocity.) Note that the intersections move to lower 
temperature with decreasing T c g. The somewhat wig- 
gly character of the i mrx curve in convection zone follows 
form the fact the the convective velocity is given by (e.g., 
Mihalas 1978) 

v = (gQH/8) 1 / 2 (V-V E )(l/H) 2 , (8) 

where g is the surface gravity, H the pressure scale 
height, Q = (din. p/dlnT)p, I is the mixing length, and 
(V — Ve) is a difference between the true temperature 
gradient, V ~ dhiT/d\nP, and the gradient of the con- 
vective elements. This is computed by considering the ef- 
ficiency of convective transport. In the standard mixing- 
length approach, the efficiencies arc computed differently 
for Rosseland optical depth smaller and larger than unity. 
The wiggly behavior corresponds to a transition between 
Rosseland optical depth larger than and smaller than 1. 
In any case, the uncertainty in the exact value of the con- 
vective velocity does not influence the non-equilibrium 
abundances of CO, CH4, and NH3 in a significant way. 

FigurefTTldepicts the abundance profiles for the models 
displayed in Fig. 1161 The non-equilibrium number frac- 
tion of CO sharply decreases mixing time with decreasing 
T c g . This is a consequence of the fact that while the in- 
tersection points of the i m ; x and the tco curves occur at 
similar temperatures around 1100 K, these intersection 
points are reached at much higher pressures for lower T e g . 
Consequently, the equilibrium number fraction of CO at 
the intersection point is lower. For ammonia, the equi- 
librium abundance profiles become flatter for deceasing 
T e fj and, therefore, the differences between the equilib- 
rium and non-equilibrium number fractions are smaller. 

In Fig. [I8l we depict the predicted flux for the three 
models displayed in Figs. [T6land[T7l The hottest model, 
at T e s = 600 K, exhibits the same features as late T 
dwarfs (see iJ3.ll) . namely the non-equilibrium flux is 
higher in the 8—14 /an region (by about 15-20%) and 
smaller in the the 4—5 /im region (on average, by about 
10%). In the cooler model at T D g = 400 K, the departures 
in the 8—14 /im region due to non-equilibrium chemistry 



of ammonia are slightly smaller than in the 600-K model. 
However, the differences in the 4—5 /im region essentially 
disappear because the non-equilibrium mixing ratio of 
CO, though much higher than the equilibrium one, is 
already quite small. Therefore, CO contributes but lit- 
tle to the opacity in this region. Finally, for the coolest 
model at T e g — 200 K, the non-equilibrium effects almost 
disappear in both the 4—5 /tm and 8—14 /mi regions. 

To diagnose the trends described above, we display in 
Fig. HU the temperature/pressure profiles for the same 
models (both non-equilibrium and equilibrium) depicted 
in Figs. [16] - [HI In Fig. [20j we display the "radia- 
tion formation temperature" as a function of wavelength 
for the two extreme models with T c g = 200 and 600 K. 
The radiation formation temperature is defined as a lo- 
cal temperature in the atmosphere where the monochro- 
matic optical depth is 2/3. This quantity indicates the 
position/layer from which a photon of a given wavelength 
effectively decouples from the atmosphere. 

4. MODEL COMPARISON WITH GLIESE 570D 
OBSERVATIONS 

In Fig. EU we display the observed spectrum of Gliese 
570D, kindly supplied to us in digital form by M. Cushing 
(Cushing et al. 2006). This is one of the latest known 
T dwarfs, and is the object studied by Saumon et al. 
(2006). We also display two theoretical models com- 
puted for T c g = 800 K and g = 10 5 cm s~ 2 , namely 
an LCE and a non-LCE model. The non-equilibrium 
model is computed with K zz = 10 cm 2 s , and for the 
"fast2" CO/CH4 reaction rate, but, as demonstrated in 
Fig. [El this predicted flux essentially coincides with the 
flux computed for all other values of K zz and CO/CH4 
reaction time prescriptions. The agreement between the 
observed and predicted spectrum for non-LCE model is 
excellent in the region between 8 and 14 /im, while the 
fit is much worse for the LCE model. We thus confirm 
the results of Saumon et al. 3 . In a future work, we will 
show comparisons between our non-equilibrium spectral 
models and observations in the near- and mid-IR of other 
late T dwarfs. 

5. DISCUSSION AND CONCLUSIONS 

In this paper, we have performed a systematic study 
of the spectral consequences of departures from chem- 
ical equilibrium in carbon and nitrogen chemistry in 
the atmospheres of L dwarfs, T dwarfs, and substellar- 
mass objects down to 200 K. Unlike previous investi- 
gations, the temperature/pressure profiles of the non- 
equilibrium models are fully consistent with the non- 
equilibrium chemistry. Such models reveal that failure 
to calculate the non-equilibrium atmospheric profiles in 
a self-consistent fashion can result in errors of as much as 
~50— 100 K in the atmospheric temperatures at a given 
pressure. Our effort has generated an extensive grid of 
non-equilibrium model spectra for effective temperatures 
from 700 K to 1800 K (in steps of 100 K), for three values 
of surface gravity, g = 10 45 , 10 5 , and 10 5 5 cm s -2 , for 
four values of the coefficient of eddy diffusion in the radia- 
tive zone K zz = 10 2 , 10 4 , 10 6 , and 10 s cm 2 s" 1 , and for 

3 Note that the predicted fluxes for both our models and those 
of Saumon et al. (2006) are too low between 7 and 8 fim, for both 
LCE and non-LCE models. The reason for this is unclear, but may 
be related to deficiencies in the extant water opacity databases. 
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three different CO/CH4 reaction time prescriptions. We 
have also provided clear and cloudy model variants. We 
find, in keeping with previous studies, that there are es- 
sentially only two regions where the effects of departures 
from chemical equilibrium significantly influence the pre- 
dicted spectrum. These arc in the M (~4— 5 /j,m) and N 
(8—14 fim) bands due to CO and NH3, respectively. The 
overabundance of CO can be very large, but translates 
into flux suppressions in the M band of at most ~40% 
between effective temperatures of 800 and 1800 K. The 
effect is largest around T e g « 1100 K. The underabun- 
dance of ammonia can be significant, but translates into 
flux enhancements in the N band of no more than ^20% 
for the T e g range from 700 to 1800 K, with the largest 
effects at the lowest values of T c g. The magnitude of 
the departure from chemical equilibrium increases with 
decreasing gravity, with increasing eddy diffusion coef- 
ficient, and with decreasing speed of the CO/CH4 re- 
action. Though these effects are modest, they lead to 
better fits with the measured T dwarf spectra (see also 
Saumon et al. 2006,2007). 

In addition, we have calculated an exploratory se- 
quence of non-equilibrium models down to 200 K, includ- 
ing the possible effects of water cloud formation. These 
models show that in the region "beyond the T dwarfs" 
(Burrows, Sudarsky, & Lunine 2003) the suppression ef- 
fect of non-equilibrium CO in the M band diminishes, 
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then disappears, with decreasing T g-. 

Hence, despite the modest suppression in the M band 
due to non-equilibrium CO chemistry for T ffS above 
~500 K, non-equilibrium effects do not compromise the 
use of the M band flux as a useful diagnostic of cool at- 
mospheres and as one of their most prominent features. 
At the higher Toil's above ~500 K, this is due in part 
to the saturation of the augmented CO absorption, but 
importantly it is due to the fact that only about half 
of the M band overlaps with the 4.7-^m CO absorption 
feature. At lower T c gs below ^400 K, non-equilibrium 
effects in the M and N bands all but disappear. Overall, 
the prominence and importance of the M band flux for 
searches for the coolest dwarfs and for brown dwarfs ( "Y" 
dwarfs ?) beyond the T dwarf range is not substantially 
affected by non-equilibrium chemistry. 
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Fig. 1. — A comparison of the equilibrium (dashed lines) and non-equilibrium (solid lines) abundances of water (black), methane 
(yellow), ammonia (red), and CO (blue), for a representative model with T c ft = 900 K, g = 10 5 ' 5 cm s — 2 , and K zz = 10 4 cm 2 s — 1 . In this 
case, water and methane remain essentially at their equilibrium abundances, while ammonia becomes underabundant at the upper layers 
(leading to a lower opacity and, thus, higher predicted flux in the 8-14 (im region for non-equilibrium models). The plots also show that 
CO exhibits larger departures from equilibrium with the slow reaction. See Fig. [5] and the text for a discussion. 



Hubeny, & Burrows 




_1 I I L 



g = 10 5 - 5 cm s-2, K = 10 4 cm 2 s" 1 ; fast2 . 



T /100 K = 8, 10, 12, 14, 16, 1! 



J I I I I I I I I I I I I I I I I I L_ 



500 



1000 



1500 



2000 



2500 



3000 



T (K) 



Fig. 2.— Similar to Fig. Q] but with a range of effective temperatures to display the trends with T e fj . The color scheme is the same 
as in Fig. [T] The association of the individual lines with effective temperatures is made by noting that in the plot the lines corresponding 
to progressively lower effective temperatures end at progressively lower local temperatures. (In other words, the longer the lines go to the 
right, the higher the effective temperature.) In all cases displayed here, water and methane remain at their equilibrium abundances. The 
degree of departure from equilibrium for ammonia decreases with increasing effective temperature. 
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Fig. 3. — Similar to Fig. [2] but showing the trend with changing surface gravity. This figure demonstrates that the non-equilibrium 
abundance of CO increases with decreasing gravity. See explanation in the text. 
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Fig. 4. — A plot of the mixing time (solid lines), the reaction time for the CO/CH4 (dashed lines), and the reaction time for the NH3/N2 
reaction (dash-dotted lines), for models with T g = 900 K, K zz = 10 4 cm 2 s^ 1 , and log 10 g (cm s — 2 ) = 5.5 (black), 5.0 (blue), and 4.5 
(red). The abrupt change of the mixing-time curves represents the boundary of the convection zone. 
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Fig. 5. — Similar to Fig. [4] but depicting the mixing times (solid lines) corresponding to (from top to bottom) the four values of the 
diffusion coefficient, K zz = 10 2 , 10 4 , 10 6 , 10 8 cm s , for models with T e g = 900 K, and g = 10 5,5 cm s . Also shown are the reaction 
times for CO/CH4 (dashed lines) chemistry, for the three different treatments of the reaction speed, "slow," "fasti," and "fast2" (the curves 
are labeled accordingly), and the reaction times for the NH3/N2 reaction (dot-dashed lines). 
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Fig. 6. — Similar to Fig. [4] but depicting the mixing times (solid lines), and the reaction times for CO/CH4 (dashed lines) for the "fast2" 
reaction rate, and the reaction time for the NH3/N2 reaction (dot-dashed line), for a set of models with g = 10 5 ' 5 cm s~ 2 , K zz = 10 4 cm 2 
s _1 , and T c ff/100 K = 9, 12, 15, 18. The reaction times for CO/CH4 are labeled accordingly. 



Non-chemical equilibrium models 



13 




Fig. 7. — A comparison of temperature-pressure profiles for equilibrium (dashed lines) and non-equilibrium (solid lines) models of 
T efI = 800, 1000, 1200, 1400, 1600, and 1800 K (from bottom to top), g is equal to 10 5 ' 5 cm s~ 2 . 
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Fig. 8. — Synoptic comparison of model atmosphere spectra computed with chemical equilibrium (black lines), and with departures 
from equilibrium (red lines), for a representative value of g = 10 5,5 cm s -2 , and K zz = 10 4 cm 2 B , for a "slow" CO/CH4 reaction, and 
for effective temperatures between 700 and 1800 K, covering most of the T and L spectral types. There are essentially two wavelength 
regions where the departure from chemical equilibrium influences the predicted spectrum: the region around 4.7 fim (due to non-equilibrium 
abundances of CO), and a wide region between 8 and 14 /im, due to non-equilibrium abundances of NH3. The CO feature has its peak for 
effective temperatures between 800 K and 1200 K; the effect generally decreases for hotter models and essentially disappears at 1800 K. 
The figure is similar to that presented by Saumon et al. (2003). 
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Fig. 9. — The same as Fig. \E\ the only difference is that the models are computed with clouds. We assume a cloud composed of a 
collection of magnesium silicate componds represented by a single species, in this case forsterite, with an extension of the cloud layer toward 
higher pressures, as suggested by BSH (E-type cloud model in their notation). We assume a modal particle size of 100 [irci. The presence 
of clouds leads to differences in the predicted spectrum, in particular for wavelengths below 2.5 )im, but otherwise the effect of departures 
from chemical equilibrium is essentially the same as for the cloudless models. 
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Fig. 10. — A comparison of predicted spectra for a model with T e ff = 900 K, K zz = 10 4 cm 2 s _1 , and three different surface gravities. 
We display the equilibrium models (black lines), and non-equilibrium models (red lines) assuming the "slow" CO/CH4 reaction rate of 
Prinn & Barshay (1977). The plot also clearly shows that the effects of departures from chemical equilibrium are more pronounced for 
lower gravities. See Fig. [T]and the text for explanations. 
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Fig. 11. — A more detailed view of the CO non-equilibrium region between 3 and 5 /im, for models with T c ff = 900 K and g = 10 5 ' 5 cm 
s — 2 . Displayed are the equilibrium model (black), as well four non-equilibrium models with K zz = 10 2 (blue), 10 4 (red), 10 6 (yellow), and 
10 s (green) cm 2 s~ 1 . All models are for the "slow" reaction rate for the CO/CH4 reaction of Prinn &; Barshay (1977). 
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Fig. 12. — Similar to Fig. 1111 but for the "fast2" CO/CH4 reaction rate of Yung et al. (1988). Departures from equilibrium are now 
smaller, as explained in the text. 
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Fig. 13. — A more detailed view of the ammonia region between 8 and 14 /im, for the same models as displayed in Fig. 1121 The curves 
for all non-equilibrium models essentially coincide, except for the highest value of K zz = 10 s cm 2 s _1 . The flux in the region between 8 
and 14 /an is relatively insensitive to the value of K zz . 
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Fig. 14. — Departures from chemical equilibrium in the 4.7-^tm CO region. The plot displays the ratio of the integrated flux in the 
region between 4.4 (im and 5 /im of the non-equilibrium and equilibrium models, as a function of effective temperature. The upper panel 
displays the dependence of the non-equilibrium ratio for g = 10 5 cm s~ 2 , for the "slow" CO/CH4 reaction (the upper set of diamonds 
connected by dashed lines) and the "fast2" CO/CH4 reaction (the lower set of stars connected by dotted lines); each set is shown for three 



values of K zz . The curves are labelled by the value of 1 



, K zz in cm 2 s" 1 . The curves for K zz = 10 6 cm 2 s" 1 for the "fast2" CO/CH4 



reaction and for K z 



s for the "slow" CO/CH4 reaction essentially coincide. The lower panel displays the dependence of the 



non-equilibrium ratio on surface gravity for two sets of models. In all models, K zz = 10 4 cm 2 s 1 . The upper set of models (diamonds 
connected by dashed lines) are for the "fast2" CO/CH4 reaction, while the lower set of models (stars connected by dotted lines) are for the 
"slow" CO/CH4 reaction. The curves are labelled by the value of log 10 g in cm s~ 2 . 
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Fig. 15. — Similar to Fig. 1141 but for the ammonia region between 9 /im and 14 (im. Since the predicted spectra are insensitive to 
the speed of the CO/CH4 reaction, and only weakly sensitive to K zz , we display only models with K zz = 10 4 , and the "fast2" CO/CH4 
reaction. The curves are labelled by the value of log 10 g in cm s — 2 . 
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Fig. 16. — A plot ol the mixing time (solid lines), the reaction time for the CO/CH4 (dashed lines), and the reaction time for the NH3/N2 
reaction (dash-dotted lines), for models with g = 10 5 cm s -2 K Z z = 10 4 cm 2 s _1 , and T c g = 600 K (black lines), 400 K (blue lines), and 
200 K (red lines). 
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Fig. 17. — A comparison of the equilibrium (dashed lines) and non-equilibriu m ( solid lines) abundances of water (black), methane 
(yellow), ammonia (red), and CO (blue), for the same models as displayed in Fig. 1161 Thin lines depict models for T c g = 600 K, thicker 
lines for for T a g = 400 K, and the thickest lines for T e g = 200 K. The curves for CO are labeled by the value of T e ff in K. 
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Fig. 19. — A comparison of temperature-pressure profiles for equilibrium (dashed lines) and non-equilibrium (solid lines) models of 
T efI = 600, 400, and 200 K (from top to bottom), g is equal to 10 5 cm s~ 2 . 




Fig. 20. — Radiation formation temperature as a function of wavelength for models with T e g = 600 K (upper curve) and 200 K (lower 
curve), g is equal to 10 5 cm s~ 2 . 
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Fig. 21. — Comparison of observed and predicted spectra of Gl 570D. Thick black histogram displays the observed flux, the red and blue 
lines are the predicted spectra for non-equilibrium (red) and equilibrium (blue) models, computed for T e ff = 800 K and g = 10 5 cm s — 2 . 
The non -equilibrium model displayed is computed with K zz = 10 4 cm 2 s _1 , and for fast CO/CH4 reaction rate, but, as demonstrated in 
Fig. 1131 this predicted flux essentially coincides with the flux computed for all other values of K zz and CO/CH4 reaction time prescriptions. 



